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Abstract. We present broad band photometry and spectroscopic observations of the host galaxy of GRB 030329. Analysis of 
the spectral emission lines shows that the host is likely a low metallicity galaxy (Z ~ 0.004). The spectral energy distribution 
(SED) constructed with the photometric points has been fitted using synthetic and observational templates. The best SED fit 
is obtained with a starburst template with an age of ~ 150 Myr and an extinction ~ 0.6. We find that the GRB 030329 
host galaxy is a subluminous galaxy (L ~ 0.0161/*) with a stellar mass of > 10* A/©. Three independent diagnostics, based 
on the restframe UV continuum, the [O II], and the Balmer emission lines, provide a consistent unextincted star formation rate 
of ~ O.6M0 yr~^, implying a high unextincted specific star formation rate (~ 34M0 yr~^ (L/L*)"^). We estimate that the 
unextincted specific star formation rate of the GRB 030329 host is higher than ~ 93.5% of the galaxies at a similar redshift. 
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1 . Introduction 

An extremely bright gamma-ray burst (GRB) was detected 
on 2003 March 29 at 11:37:14.67 UT (GRB 030329; 
Vanderspek et al. 2003 ) by the French Gamma Ray Telescope 
(FREGATE), the Wide Field X-Ray Monitor (WXM), and 
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* Based on observations made with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astroffsica de Canarias. Based 
on data taken at the 2.2-m and 3.5-m telescopes of the Centro 
Astronomico Hispano Aleman de Calar Alto, operated by the Max 
Planck institute of Heidelberg and Centro Superior de Investigaciones 
Cientfficas. The spectral observations were obtained at the European 
Southern Observatory, Cerro Paranal (Chile), under the Director's 
Discretionary Time programme 27I.D-5006(A). 



the Soft X-Ray Camera (SXC) instruments aboard the 
High Energy Transient Explorer II (HETE-2) spacecraft 
(Ricker et al. 2003). With a duration of approximately 30 s 
in the 30^00 keV energy range, GRB 030329 falls into the 
"long-duration" category of GRBs (Kouveliotou et al. ll993> . 

Before the HETE-2 era, two main scenarios were favoured 
to explain the nature of long GRBs. One of them suggested that 
GRBs occur when two compact objects (such as black holes or 
neutron stars) merge (Eichler et al. 1989 Narayan et al. ll992> . 
This would imply that GRBs occur after massive star forma- 
tion has ended in a galaxy, since the time-scale for binary ob- 
jects to merge is long (> 1 Gyr). A second scenario, the so- 
called collapsar scenario, suggested that GRBs are associated 
with the core-collapse of massive stars (supernovae or hyper- 
novae: Woosley 1993 Paczyriski 1998 MacFadyen & Woosley 
II999I 1. This scenario predicts a supernova (SN) accompanying 
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the GRB, which thus must be coincident with the epoch of star 
formation in the host. 

The studies carried out in 2003 for GRB 030329 (Hjorth 
et al.'200T, Stanek et al. l2003l Sokolov et al. 120031 Kawabata 
et al. 2003 1 strongly confirmed previous evidences (Bloom et 
al. [T999l Castro-Tirado & Gorosabel [T999t that long GRBs are 
physically linked to supernova explosions. The origin of the 
short GRBs is still unknown, so here we will discuss only the 
hosts of long GRBs. Radio, optical and/or infrared afterglows 
have been observed for ^^70 long GRBs, and the majority of 
these coincide with star forming galaxies. 

The pioneering statistical work applying the spectral energy 
distribution (SED) fitting technique to GRB host galaxies was 
published by Sokolov et al. ( 2001 1, who argued that GRB hosts 
correspond to starburst galaxies. A number of other comple- 
mentary studies have been carried out in the past, studying the 
properties of different GRB host samples (Chary et al. 2002,, 
Berger et al. l2003l Le Floc'h et al. l20U3l Tanvir et al. l2004t or 
centred on a particular host galaxy (Djorgovski et al. 1200 lal 
Djorgovski et al. 2003 Hjorth et al. 2002 1. Recently, a sta- 
tistical analysis by Christensen et al. ( 2004al hereafter C04A) 
based on the SED of a sample of 10 GRB host galaxies in- 
dicated the similarity of GRB hosts and the young starburst 
galaxies present in the Hubble Deep Field (HDF; Williams et 
al.EOOO^. 

A good knowledge of the host SED makes it possible to 
derive the key parameters for determining the conditions that 
favour a GRB event, e.g., the star formation rate (SFR), the 
dominant stellar population age, the stellar mass and the extinc- 
tion. GRB host galaxies tend to be faint (R > 23; Djorgovski 
et al. 1200 IM. and thus, spectroscopic studies of the SED re- 
quire extensive use of large telescopes. An economic alterna- 
tive that allows us to extract information on the properties of 
the host galaxies is based on multicolour broad band imag- 
ing. The main inconvenience of fitting templates to photomet- 
ric points comes from the insensitivity of the SED solution to 
the assumed metallicity (Christensen et al. '2004b^, which is 
considered as a secondary parameter (Bolzonella et al. 2000 1. 
In practice, in order to determine reliable metallicities for the 
GRB hosts, spectroscopic observations are therefore required. 

The optical afterglow (OA) coordinates (aj20oo — 
10''44™59!95, (5j20oo = +21°31'17'/8_with an uncertainty of 
0.5 arcsec in each axis; Price et al. '2003^ allowed an exten- 
sive follow up from the Northern hemisphere, although it was 
also visible from the main observatories of the Southern hemi- 
sphere. The redshift of GRB 030329 was determined to be 
z = 0.168 (Greiner et al. 2003 1 from high-resolution spec- 
troscopy with the Very Large Telescope (VLT). This makes 
GRB 030329 the third nearest burst overall (GRB 980425 
is the nearest at z = 0.0085, Galama et al. i998, and 
GRB 031203 had z = 0.1055, Prochaska et al. l2054t . As the 
OA faded, prominent spectral features revealed a Type Ic super- 
nova (named SN 2003dh) peaking 10-15 days after the GRB 
event (Hjorth et al. 2003, Stanek et al. 2003 ). 

In this paper we present a comprehensive multicolour study 
of the GRB 030329 host galaxy, similar to those performed 
for the GRB 000 210 (Go rosabel et al. 2003al, GRB 000418 
(Gorosabel et al. I2003b> and GRB 990712 (Christensen et 



al. 2004bl host galaxies. The aim of the analysis is to deter- 
mine the properties of the stellar populations dominating the 
optical/near-infrared (NIR) light from the host galaxy and the 
amount of extinction due to dust in the interstellar medium 
(ISM) of the host. 

Given the faintness of GRB host galaxies in general, de- 
tailed observations of the low redshift hosts, as in the present 
study, might be used to construct a library of spectral templates, 
applicable (under several assumptions) to the hosts placed at 
higher redshifts. Furthermore, the case of GRB 030329 is par- 
ticularly remarkable thanks to its association with SN 2003dh. 
Both arguments justify a detailed analysis of this host galaxy. 

Throughout this paper, we assume a cosmology where 
Ha = 0.7, = 0.3 and Hq = 65 km s'^ Mpc~^. Under 
these assumptions, the luminosity distance of GRB 030329 is 
di = 872.5 Mpc and the look-back time is 2.25 Gyr (15.5 % of 
the present Universe age). 

Sect. IJldescribes the photometric and spectroscopic obser- 
vations carried out for the GRB 030329 host galaxy. Our main 
results are reported in Sect.|3] In Sect.l^we discuss the implica- 
tions of our work and Sect.|5]summarises the final conclusions. 

2. Observations and Analysis 

2.1. Imaging 

We have used a number of optical/NIR facilities in order to 
compile a well sampled SED (see Table [Q. The UBVR- 
bands were observed with the 2.56-m Nordic Optical Telescope 
(NOT) equipped with MOSCA, a 2x2 mosaic 2048x2048 
pixel CCD, covering a field of view (FOV) of 7.'7 x 7.'7. These 
observations were carried out in a 2x2 binning mode, yielding 
a pixel scale of O'.'22/pix. 

Additional optical observations were obtained with the 
BUSCA camera (Reif et al.j^OOOl at the 2.2-m telescope lo- 
cated at CAHA'. BUSCA allows simultaneous imaging in four 
broad optical bands. The four channels (hereafter named CI, 
C2, C3 and C4) are not standard filters and have to be cali- 
brated by observing spectro-photometric standard stars (char- 
acteristics of the four BUSCA channels are described in Table 

CI resembles the Johnson [/-band, C2 is a transition be- 
tween Johnson's B and V filters. C3 is very close to the Rc- 
band and C4 is basically identical to the Ic filter Figure [0 
shows a coloured composite image based on three out of four 
BUSCA channels, namely CI, C2, and C4. BUSCA consists 
of four 4Kx4K pixel CCDs plus 3 dichroics. The observations 
were performed in the 2x2 binning mode. The FOV covered is 
12f0 X 12f0, and the resulting pixel scale was O'.'35/pix. 

Finally, the JHK' observations were carried out with 
the 3.5-m CAHA telescope equipped with Omega-Prime, a 
1024x1024 pixel HgCdTe CCD, which provides a FOV of 
6f8 X 6.'8 and a pixel scale of G'.'39/pix. 

Reduction was performed following standard procedures 
running under IRAF^. In Table^we provide the observing log 
for our optical and NIR imaging observations. 

' Centre Astronomico Hispano Aleman 

" IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
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Telescope 


Filter 


Date UT 


^cxp 


Seeing 


Limiting magnitude 


(+lnstrument) 






(s) 


(arcsec) 


(3a) 


NOT (+MOSCA) 


U 


24.03-24.08/03/04 


5x900 


1.0 


25.4 


NOT (+MOSCA) 


B 


24.85-24.89/03/04 


3x600 


1.1 


25.8 


NOT (+MOSCA) 


V 


24.86-24.87/03/04 


1x600 


1.0 


24.4 


NOT (+MOSCA) 


R 


24.87-24.88/03/04 


1x600 


0.9 


24.2 


2.2-m (+BUSCA) 


CI 


20.94-21.01/03/04 


19x600 


1.8t 


25.6* 






21.97-22.12/03/04 


14x600 










22.95-23.06/03/04 


12x600 






2.2-m (+BUSCA) 


C2 


20.94-21.01/03/04 


19x600 


1.7t 


25.1* 






21.97-22.12/03/04 


14x600 










22.95-23.06/03/04 


12x600 






2.2-m (+BUSCA) 


C3 


20.94-21.01/03/04 


19x600 


1.7t 


25.2* 






21.97-22.12/03/04 


14x600 










22.95-23.06/03/04 


12x600 






2.2-m (+BUSCA) 


C4 


20.94-21.01/03/04 


19x600 


1.6t 


24.1* 






21.97-22.12/03/04 


14x600 










22.95-23.06/03/04 


12x600 






3.5-m (+Omega-Prime) 


J 


05.07-05.19/01/04 


113x60 


0.7 


24.1 


3.5-m (+Omega-Prime) 


H 


06.07-06.18/01/04 


109x60 


0.8 


22.9 


3.5-m (+Omega-Prime) 


K' 


07.06-07.20/01/04 


99x60 


0.8 


19.7 



f Mean seeing of the co-added image. 

★ Limiting magnitude obtained co-adding the three epochs. 



Table 2. Optical/NIR magnitudes of the host. The fourth and fifth columns show the measured magnitudes in the Vega and AB 
systems (not corrected for Galactic reddening). 



Filter name 


Effective 


Bandpass 


Magnitude 


Magnitude 


(Instrument) 


wavelength (A) 


width (A) 


(Vega) 


(AB) 


U (MOSCA) 


3638.3 


396.3 


22.72±0.10 


23.45±0.10 


CI (BUSCA) 


3654.5 


1081.6 


23.03±0.03 


23.51±0.03 


B (MOSCA) 


4358.7 


680.5 


23.33±0.07 


23.26±0.07 


C2 (BUSCA) 


5007.4 


706.2 


22.92±0.05 


22.87±0.05 


V (MOSCA) 


5413.4 


636.4 


22.83±0.10 


22.84±0.10 


R (MOSCA) 


6422.1 


1072.5 


22.66±0.04 


22.86±0.04 


C3 (BUSCA) 


6455.8 


1147.0 


22.61±0.04 


22.81±0.04 


C4 (BUSCA) 


7999.7 


780.7 


22.10±0.04 


22.55±0.04 


J (Omega-Prime) 


12861.5 


1713.9 


21.45±0.16 


22.42±0.16 


H (Omega-Prime) 


16491.6 


1814.2 


21.15±0.24 


22.55±0.24 


A"' (Omega-Prime) 


21634.9 


1636.5 


>19.70t 


> 21.57t 


f 3 (7 upper limit 



The optical/NIR imaging was carried out at least 282 days 
after the GRB, when the OAh-SN contamination was negligible. 
The host galaxy appears as a star-like object in all our images, 
showing NIR and optical profiles consistent with those of field 
stars. This is consistent with the HST observations by Fruchter 
et al. J2003> . who reported that the host is a dwarf galaxy with 
a brightness of ^ = 22.7 ± 0.3 (consistent with our T^-band 
magnitude, see Table |2}. According to the angular scales of the 
OA-host system discussed in Fruchter et al. ( 2003 ), a host com- 
pact appearance is expected for seeing conditions above ^ 0'.'5, 

for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



in agreement with our observations (seeing > Cf.'l; see Table0. 
Therefore, considering that the relative photometry is indepen- 
dent on the aperture radius, the host magnitudes given in Table 
12] are based on circular aperture measurements (PHOT running 
under IRAF) with no aperture corrections and an aperture ra- 
dius of 1.5 times the full width at half maximum (FWHM). 

The U BVR-hand calibration is based on Henden J2003> . 
and the JHK' calibration on Guziy et al. (I2005> . For the 
BUSCA data, the flux calibration was caiTied out by observ- 
ing the spectro-photometric standard star GD153 (Bohlin et al. 
1995 1 at the same airmass to that of the GRB host on a photo- 
metric night. This allowed us to establish secondary field stan- 
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Fig. 1. The field of the GRB 030329 host galaxy. The pic- 
ture is a coloured composite image created by combining three 
BUSCA images taken in the CI, C3, and C4 bands. The FOV 
is 2'.0 X 2'.0. The location of the host galaxy is indicated. As 
can be seen from the picture, the host shows a bluish colour 
North is upwards and East is leftwards. 



dard stars in the AB system. As can be seen in Fig.|3] the agree- 
ment of the BUSCA calibration with the NOT optical magni- 
tudes is satisfactory. 

2.2. Spectroscopy 

Spectroscopy of the host of GRB 030329 was obtained on 
2003 June 19 with the Focal Reducer and Low Dispersion 
Spectrograph (FORS2) instrument^ attached to the 8.2-m UT4 
of the European Southern Observatory Very Large Telescope. 
The FORS2 instrument is equipped with two 2048 x 4096 pixel 
MIT CCDs having 15 /im pixels. We used the standard resolu- 
tion mode that gives a FOV of 6'.8 x 6.'8, and with 2x2 binning 
on the chip we obtained a pixel size on the sky of 0'.'25. 

The observations were conducted with the 300V grism and 
an order sorting filter GG375, which efficiently covers the 
wavelength range from ~ 3800-8800 A. A 1'.'3 wide slit was 
used and the seeing during the observations was typically be- 
low 0'.'6, yielding a spectral resolution of ^ 10.5 A. Given the 
small size of the object, the entire flux of the host galaxy should 
thus be included in the slit. A position angle of 123.6 degrees 
was used to align also a nearby star on the slit (see Hjorth 
et al. I2003> . Three spectra of 900 s each were obtained, with 
small offsets on the slit between the individual exposures. The 
airmass was always quite substantial since these observations 



were obtained from the Southern hemisphere, and varied from 
1.8 to 2.3. The final combined spectrum is shown in Fig.|2] 

The spectra were reduced in a standard way, including 
bias subtraction, flat fielding, and wavelength calibration us- 
ing spectra of a Helium- Argon lamp. Flux calibration was done 
relative to the spectrophotometric standard star Feige67 (Oke 
I990i. The absolute flux calibration of the combined spectrum 
was obtained relative to our broad band photometry. 

We note that our host galaxy spectroscopic observations 
were carried out ^ 82 days after the GRB, when the contri- 
bution of the OAh-SN and the host were comparable (Guziy et 
al. l2005> . In contrast, the host properties derived by Matheson 
et al. (2003'l are based on spectra acquired only 6 — 10 days af- 
ter the GRB, when the OAh-SN was still 50 - 25 times brighter 
than the host. This intense overimposed radiation field might 
affect the measurement of the host emission lines fluxes (es- 
pecially the faint ones), so we consider our spectrum a cleaner 
probe of the host galaxy physical properties. 

Before using the emission lines as indicators of several 
host galaxy properties, they were corrected for Galactic fore- 
ground reddening in the direction of the GRB (E{B ~ V) = 
0.025; Schlegel et al. lI998> considering a MW extinction law 
(Cardellietal. .l989i l. 

2.3. Photometry: Constructing the optical/NIR SED 

In order to derive the corresponding effective wavelengths and 
AB magnitudes we convolved each filter transmission curve 
plus the corresponding CCD efficiency curve (see Table|2} with 
the Vega spectrum (a Lyrae, m = in all bands by definition; 
Bruzual & Chariot 1993 1. The fit of the SED has been mainly 
obtained using synthetic templates, but checks have also been 
performed using observational templates (Kinney et al. 119961 
hereafter K96). 

The synthetic SED analysis is based on the code HyperZ"* 
(Bolzonella et al. |2000l. Eight synthetic spectral types were 
used representing Starburst galaxies (Stb), Ellipticals (E), 
Lenticulars (SO), Spirals (Sa, Sb, Sc and Sd) and Irregular 
galaxies (Im). Each spectral type has associated a SFR(i) tem- 
poral history, and hence a characteristic time scale (SFR cx 
exp (— t/r); where t can range from for Stb, to cxd for Im 
galaxies). For the generation of the synthetic templates, two 
initial mass functions (IMFs) have been considered; Salpeter 
(|1955; S55), and Chabrier ( l2U03l CH03). Four extinction laws 
have been taken into account: Calzetti et al. (2000; suitable 
for Stbs), Seaton iWM for the Milky Way, MW), Fitzpatiick 
( 119861 for flie Large Magellanic Cloud, LMC), and Prevot et 
al. ( fT984l for the Small Mageflanic Cloud, SMC). The con- 
struction of the HyperZ templates has been performed using 
the GALAXEV public code (Bruzual & Chariot 2003), using 
as input the metallicity (Z, see Sect. 13 .4t derived from our VLT 
spectroscopic observations. 

In addition to the aforementioned evolutionary templates, 
the observational templates from K96 were considered. These 
extra spectral templates give complementary information on 
the galaxy type and extinction. The observed K96 templates 
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can be grouped into seven sets: Bulge (B), E, SO, Sa, Sb, Sc 
and Stb galaxies. Additionally, the Stb K96 templates are sub- 
divided into six classes depending on their extinction by in- 
creasing E{B — V) in steps of 0.1 (from Stbl to Stb6; see 
more details in K96). The extinction for the B, E, SO, Sa, Sb, 
and Sc templates is a free parameter 

For the construction of the photometric SED all the magni- 
tudes were corrected for Galactic foreground reddening along 
the GRB line of sight (Schlegel et al. ll998> assuming a typical 
MW extinction law (Cardelli et al. ll989> . 

3. Results 

The VLT spectrum allows us to estimate the host galaxy extinc- 
tion (^v), the metallicity (Z) and the SFR, based on the oxygen 
and Balmer emission lines. Table|3ldisplays the list of emission 
lines detected in our VLT spectrum. The flux of the lines, with 
the exception of H6, are consistent within 1. 5 a with the fluxes 
previously reported for the GRB 030329 host galaxy (Hjorth et 
al. l2003ll . 

Alternatively, the SED constructed from the optical/NIR 
photometric points gives an independent estimate of the ex- 
tinction and the SFR. The SED also provides us with informa- 
tion on the stellar population age, the stellar mass, and on the 
host galaxy absolute luminosity. For self-consistency, we based 
our SED fitting analysis on templates constructed adopting the 
metallicity derived from our spectroscopic study (Z — 0.004). 
Table |2] shows the optical/NIR magnitudes given in the Vega 
and AB systems, not corrected for Galactic or host reddening. 
The listed magnitude errors account for both the zero point and 
the statistical errors. Figure|3lshows our photometric points and 
the SED solution obtained for a SMC extinction law, a S55 IMF 
and a metallicity of Z ^ 0.004. Considering all the possible 
extinction laws (4) and IMFs (2), we end up with 8 possible 
templates for each galaxy type. 

3. 1. Photometric redshift 

Given that the spectroscopic redshift is known, we only use the 
photometric redshift to check our fitted SEDs for internal con- 
sistency. The spectroscopic redshift falls within the 68% un- 
certainty region of the photometric redshift in 6 out of the 8 
IMF/extinction law combinations. 

A rough estimate based on the K96 templates yields a pho- 
tometric redshift around z ^ 0.11, reasonably close to the 
spectroscopic redshift. We consider the inferred photometric 
redshifts to be basically consistent with the spectroscopic red- 
shift, and hereafter we fix the redshift of the SEDs to z = 
0.168. 

3.2. Absolute magnitude 

Subluminous galaxies are defined as having luminosity values 
below the knee of the luminosity function given by L* (or the 
equivalent Mg). Slightly different values for Mg have been re- 
ported depending on the galaxy colour and type (e.g., Schechter 
119761 Lilly et al. il995i) . For simplicity and in order to perform 



a comparative study with the C04A results, we have assumed 
the same value of A/^ = —21 as adopted by C04A. 

According to the synthetic templates, we estimate a mean 
S-band absolute magnitude of Mb = —16.5, similar to that of 
the SMC. This value is virtually independent on the assumed 
extinction law, IMF and metallicity. The Mb = —16.5 value 
corresponds to a luminosity of L ^ 0.016-L*. We conclude that 
the GRB 030329 host is clearly a subluminous galaxy. This is 
consistent with the host luminosity inferred by Fruchter et al. 

3.3. Estimating 

Among the 8 possible extinction law and IMF combinations, 
we find for 50% of the cases an extinction value within the 
0.6 < Ay < 0.7 bin. For these four combinations no clear 
preferential extinction law or IMF is found. The rest of the in- 
ferred Av values are equally distributed in the 0.7 < A^ < 0.8 
and Av < 0.6 regions. The mean extinction value derived for 
the synthetic templates is A^ — 0.6. 

We have also used the empirical templates by K96 to verify 
the extinction value obtained with the synthetic templates. The 
best fit is then obtained with the blue Stbl and Stb2 SEDs (see 
Table 13 X^/d-o.f ^ 3), which are defined to have the low- 
est extinctions among the K96 Stb sample (0 < E{B — V) < 
0.21). Stb templates with larger E{B — V) values clearly yield 
worse results. To translate the value of E{B — V) to A^, we 
have adopted a value of Ay/E{B - V) = 3.1. Thus, the Stbl 
and Stb2 templates imply a value of < Av < 0.65, in agree- 
ment with the mean A^ = 0.6 inferred based on the synthetic 
templates. The assumed Ay/E{B — V) fraction is indicative 
(the actual Ay/E{B — V) value depends on the unknown host 
extinction law) and only aimed to check the A^ value inferred 
with the synthetic templates. 

Figure |3] shows the similarities of the Stbl empirical 
template and the solution obtained by an evolutionary SED. 
Leaving E{B — V^) as a free parameter for the non Stb tem- 
plates, we obtained the E{B — V) values reported in column 
3 of Table These fits are not satisfactory. Higher values 
of E{B — V) for these templates would make the spectrum 
even redder and hence provide even larger /d.o.f. values. 
The exception is the Sc template which presents a value of 
/d.o.f. = 3.47, just slightly larger than for the Stbl and 
Stb2 templates. The E{B — V) inferred for the Sc template is 
actually consistent with the one derived for the empirical Stbl 
and Stb2 templates, as well as with that derived for the syn- 
thetic SEDs. 

In principle, the relative fluxes of the Balmer lines can be 
used as an indicator of the intrinsic host galaxy A^ (Osterbrock 
1989). The Ha emission line is affected by a nearby telluric 
band, so we have used the rest of the Balmer lines present in 
our spectrum (li/3, H7, HS) to estimate A^. A major difficulty 
in using the faint Balmer lines comes from the potential con- 
tribution of an underlying evolved stellar population. This can 
introduce absorption features with equivalent widths (EWs) up 
to a few A (Bruzual & Chariot 2003 ), which are comparable to 
our measured emission fluxes in H7 and HS. We have averaged 




4000 4500 5000 5500 6000 6500 7000 7500 



Observed wavelength (A) 

Fig. 2. The GRB 030329 host galaxy VLT spectrum. The flux is corrected for Galactic reddening {E{B — V) = 0.025) consid- 
ering a MW extinction law (Cardelli et al. 1989). The residuals resulting from the telluric line subtractions are indicated by an S. 
Note the telluric line affecting the left wing of Ha. The spectrum was acquired ^ 82 days after the GRB, when the contribution 
to the continuum by the OA+SN was still significant. 



the H7 and US absorption EWs of the 8 galaxy templates fit- 
ted to our photometric points and corrected the emission lines 
fluxes to estimate A^. The mean H7 and HS restframe absorp- 
tion EWs were 8±2 A and 7±3 A, respectively. The extinction 
values inferred are = 0.33 ± 0.19 and A^ ^ 0.86 ± 0.22, 
based on H7 and HS, respectively. The weighted mean of the 
two extinction values is A^ — 0.56 ± 0.14, consistent with 
the j4v — 0.6 value previously inferred based on the synthetic 
and empirical galaxy templates. We note that the A^ values 
determined through H7, and especially from HS, are signifi- 
cantly sensitive to the correction for the underlying population. 
Therefore, although the agreement is very good, we will as- 
sume the extinction determined with the synthetic SED tem- 
plates and use A^ — 0.6 throughout the paper. 



Table 3. Emission lines detected in the GRB 030329 host 
galaxy VLT spectrum. The fluxes are corrected for Galactic ex- 
tinction. 



Element Flux Observed Redshift 



(10 erg cm ^) wavelength (A) 



[On] 


20.40 ± 0.39 


4354.08 


0.1681 


[Nelll] 


4.04 ±0.28 


4520.15 


0.1682 


H(5 


2.33 ± 0.20 


4792.96 


0.1684 


H7 


4.82 ± 0.20 


5071.08 


0.1683 


H/3 


11.70 ±0.36 


5679.70 


0.1683 


[Oiii] 


14.76 ±0.38 


5793.80 


0.1684 


[Oiii] 


42.21 ±0.31 


5849.79 


0.1684 


Ha 


32.40 ± 0.47 


7667.55 


0.1683 


[Sii] 


3.27 ±0.32 


7861.65 


0.1684 



3.4. Estimate of the metallicity 

The detections of [Oil], [Olll] and H(3 in emission allow us 
to estimate the host metallicity based on the R23 iterative tech- 
nique (e.g., Kewley & Dopita "2002"). Correcting the emission 
line fluxes for A^ ~ 0.6 and assuming a SMC extinction law, 
we obtain a metaflicity of Z = 0.004 (Z = 0.20^©, where 
Zq — 0.02) and an ionization parameter q = 6.84 x 10^ cm 
s^^. The inferred metallicity is well below Zq, implying that 
the GRB 030329 host is a low metallicity galaxy. 

The inferred metallicity is basically insensitive to the as- 
sumed extinction law. Typical MW and LMC extinction laws 



(Pei ll992l yield metallicity values only ~ 3% appart from the 
metallicity derived based on the SMC. 

Hjorth et al. (i2003i reported a rough metallicity value of 
[0/H]= -1.0 (equivalent to Z = O.IOZ©), not correcting for 
Galactic reddening, intrinsic host galaxy reddening and under- 
lying population. If the above mentioned corrections are not 
applied to our measurements, we obtain a metallicity value of 
[0/H]= -0.8 (Z = 0A6Zq). Given that the typical metalHc- 
ity uncertainty of the R23 technique is ±0.2 dex (Kewley & 
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Dopita"2002^, we conclude that our metallicity estimate is in 
agreement with the value inferred by Hjorth et al. (2003 1. 

The i?23 diagnostic is bivaluated (lower and upper branch), 
providing two possible metallicity values. In this study we as- 
sume the metallicity inferred based on the R23 lower branch. 
According to Sollerman et al. J2005t the R23 upper branch 
yields a metallicity value of [0/H]=: —0.1. This metallicity de- 
generacy can be broken by considering other emission lines 
such as [N 11] which was not detected in our VLT spectrum. 
As noted by Sollerman et al. ( I2005> . the [N 11] flux upper limit 
derived from the VLT spectrum is not constraining enough to 
resolve the R23 degeneracy. 

However, the R23 upper branch predicts a rather bright 
galaxy (Mb ~ —19; Lee et al. .2003il . while the lower branch 
metallicity is more consistent with the Mb magnitude of the 
host (Mb ^ —16 is expected for Z = O.2OZ0), as pointed 
out by Sollerman et al. (12005 1. Therefore, the lower branch so- 
lution is more likely than the upper one, but a host metallicity 
close to solar abundance can not be completely excluded. 

3.5. Estimate of the SFR 

Following the standard Ha and [Oil] diagnostic techniques 
(Kennicutt 1998, hereafter K98) the SFR can be estimated. 
However, given the contamination that Ha suffers from telluric 
lines, we have used the SFR given by the [O ll] diagnostic. This 
gives a SFR=0.26 Mq yr^^ (corrected for Galactic extinction 
but not for host galaxy extinction), in agreement with the SFR 
derived by Hjorth et al. flOOf. SFR=0.20±0.07 Mq yr" ^ once 
it is corrected for Galactic extinction). For a Pei J1992> SMC- 
like extinction law, and adopting = 0.6, we obtain an un- 
extincted SFR of 0.6 Mq yr"^ 

A verification of the SFR inferred from the [O ll] line can be 
carried out estimating three independent SFRs through H/3, HS, 
and H7 assuming a typical case B recombination and taking 
account the extinction previously derived. The inferred three 
SFRs show a mean weighted value of 0.63 Mq yr~^ and a 
dispersion of 0.15 Mq yr^^. This is fully consistent with the 
SFR derived based on [O ll]. 

The photometric SEDs bring an additional diagnostic to es- 
timate the SFR. The method for deriving the SFR from the UV 
continuum flux is one of several diagnostic methods used in the 
literature to measure SFRs in galaxies (see K98 for a compre- 
hensive review). Obviously, if there is dust-enshrouded star for- 
mation, then this UV-based method will only provide a lower 
limit to the actual SFR. At z = 0.168 the 2800 A region is red- 
shifted to ^^3270 A, so it is included in the sensitivity region of 
the J7-band and the BUSCA CI band. We have derived the av- 
erage value of the fluxes in both filters (F = 1.5 x lO^'^^ 

2800A 

erg s~^ Hz~^ cm~^) and determined the UV luminosity at rest- 
frame 2800 A to be Luv ~ 1.2 xlO^'^ erg s^^Hz^i. Following 
K98 this corresponds to a SFR=0.17 Mq yr"^ Given that the 
rest frame luminosity of the host is L = 0.016L*, we ob- 
tain a specific star formation (SSFR) rate of SSFR=10.6 AIq 
yr-i(L/L*)-i. If we coiTect the previous UV SFR and SSFR 
estimates for the reddening effect in the host galaxy (A^ ^ 0.6, 
assuming a SMC extinction law), we obtain unextincted SFR 



- 0.54 Mq yr-i and SSFR- 34 Mq yr'^{L/L*)-^. Thus, 
the unextincted SFR derived from the UV continuum is in 
agreement with the unextincted SFR inferred based on the 
[On] and the Balmer lines (SFR-- O.6Af0 yr~^). Given that 
C04A used the UV diagnostic technique for their SSFR esti- 
mates, for comparative purposes we adopted the SSFR derived 
from the UV continuum (- 34 Mq yx-'^{L/ L*)~^). 

3.6. A potential galaxy association at z 0.17? 

A GRB 030329 field galaxy located at a.72000 = 
10'H5'"01.46", (Sj2ooo = +21°29'30'.'l (with an internal 
error of 0'.'7 per coordinate) was accidentally placed on the 
FORS2 slit. The galaxy showed prominent emission lines 
similar to the host galaxy, and turned out to be at a redshift of 
z = 0.1710 ± 0.0003. The similarity to the GRB host galaxy 
redshift, naturally opens the question on a potential clustering 
or association at z — 0.17. The nearby galaxy is located 192.4 
arcsec from the host, implying a distance of — 600 kpc at the 
host galaxy redshift. The radii of galaxy clusters is of the order 
of a few Mpc, so the measured angular distance to the host 
galaxy is compatible with a physical association at z — 0.17. 
On the other hand, the velocity difference is — 685 ± 104 km 
s^^, in agreement with the velocity dispersions measured in 
z < 0.15 galaxy clusters (Desai et al. l2004> . 

A spectroscopic study of a significant galaxy sample in the 
GRB field would be necessary to investigate whether the host 
is physically linked to other angularly close galaxies. If such 
an association is confirmed, the case of the GRB 030329 host 
would be opposite to the preferential tendency of GRB hosts 
to reside in low galaxy density environments (Bornancini et al. 

IMmIi . 

3. 7. Galaxy type 

When using the synthetic templates for the fit, we found that 
the favoured galaxy type is not very dependent on the assumed 
IMF and extinction law. Among the 8 IMF/extinction law com- 
binations, a starburst galaxy is obtained in 50% of the cases. 
The rest of cases are distributed in Sa (25%), SO (12.5%) and 
Sb (12.5%). Thus the synthetic SED fits favour a Stb type. 

This result is supported by the fits obtained with the empir- 
ical templates by K96 (Table 1^. In general the starburst tem- 
plates give lower values for x^/d-o.f. when fitting the SED. 
The only exception is for the Sc type, /d-o.f.=3.47 , whose 
value is quite close to the best ones obtained for the Stbl and 
Stb2 templates. 

3.8. Dominant stellar population age 

The inferred mean stellar population age for the 8 
IMF/extinction law combinations is 150±80 Myr In principle 
progenitors older than —50 Myr are not obvious to accommo- 
date within the collapsar scenario, since the age of an ~ 8Mq 
star (the minimum SN progenitor mass) when it explodes as a 
core-collapse SN is —50 Myr (see for instance Portinari et al. 
fT99l . 
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Fig. 3. The optical/NIR SED of the GRB 030329 host galaxy. 
The circles show the UBVR, CI, C2, C3, C4 and JH-hand 
detections. The filled triangle represents the K' upper limit. 
The abscissa is represented on a logarithmic scale in order to 
enlarge the region around the 4000 A break. The dashed line 
shows the fit obtained assuming a metallicity of Z — 0.004, 
a S55 type IMF and a Prevot et al. |1984i extinction law. Note 
the similarity with the Stbl K96 empirical template overplotted 
with a continuous line. 



Table 4. Main results obtained based on the empirical tem- 
plates by Kinney et al. f 19961. The templates have been or- 
dered according to their colour, bluer on the top (Stbl) and 
redder at the bottom (B). There is a clear correlation between 
the goodness of the fit and the colour of the template: the bluer 
the colour the lower /d.o.f.. 



Galaxy 
Template 


xVd.o.f. 
(d.o.f. = 9) 


E{B - V) 


Stbl 


3.06 


0.05 ±0.051" 


Stb2 


3.01 


0.15 ±0.051" 


Stb3 


14.86 


0.30 ±0.05"!" 


Stb4 


24.48 


0.45 ±0.05"!" 


Stb5 


36.77 


0.55 ±0.05"!" 


Stb6 


47.57 


0.65 ±0.05"!" 


Sc 


3.47 


0.18 


Sb 


133.38 





Sa 


153.72 





SO 


176.87 





E 


182.18 





B 


180.96 






f Extinction fixed by the template definition. 



However, an ideal starburst where all the star formation oc- 
curs simultaneously, and hence all the stars have exactly the 
same age, is not realistic. Thus, in a real starburst, like the 
GRB 030329 host, although the majority of the star forma- 
tion occurred ^ 150 Myr ago, the star formation can not stop 



sharply, but has a stellar age dispersion which spreads down to 
150 Myr. Hence, an age of ^ 150 Myr should be considered an 
upper limit of the progenitor age. 

An independent check of the dominant stellar population 
age can be performed comparing the restframe Hf3 equiva- 
lent width derived from the VLT spectrum (E'W^SO A, cor- 
rected for the underlying population) to the one predicted by 
Starburst99 (SB99, Leitherer et al. 1999 1. According to SB99, 
for a metallicity of Z ^ 0.004 our H/3 E'W corresponds to ages 
ranging from ~ 10 Myr (for the case of an ideal instantaneous 
SFR, T = 0) to - 300 Myr (for a constant SFR, r = oo). Our 
age of 150 ± 80 Myr is consistent with this age range, sup- 
porting that the star formation might be brief, but not perfectly 
instantaneous (r = 0). 

3.9. Host galaxy stellar mass 



Given that the Bruzual & Chariot (I2003> templates are nor- 
malized to one Solar mass, the normalization factor used to 
scale the templates to our photometric fluxes allowed us to es- 
timate the host galaxy stellar mass (M^ost)- We caution how- 
ever that the optical photometric points are likely dominated 
by young massive stars, so this estimator tends to be insensi- 
tive to the redder low-luminosity (but very abundant) low-mass 
underlying stellar population. This effect is especially relevant 
in galaxies, like the GRB 030329 host, where the optical light 
is dominated by a bright starburst episode(s). So, strictly speak- 
ing, the derived mass should be considered a lower limit to the 
total stellar mass present in our host. 

Assuming the S55 and CH03 IMFs and not accounting for 
the effect of the internal extinction, we derived mass values 
of Mhost = 9.1 ± 1.3 X IOHIq and Afhost = 6.0 ± 0.9 x 
IO^Mq, respectively. Given that an extinction value of = 
0.6 implies a scale up of the intrinsic optical flux by a factor of 
^ 2, we adopt 10* Af© as a conservative lower limit to the total 
stellar mass. 

An independent rough estimation of this value can be ob- 
tained through the statistical mass-to-luminosity ratios {M/L) 
reported for the galaxies of the Sloan Digitized Sky Survey 
(SDSS; York et al.'2000'). Adopting K = K', for the i^T-band 
absolute magnitude {AIk = —19.1, basically not affected by 
the host galaxy reddening if ^ 0.6 is assumed) and the 
B — R restframe colour of our host, M/L ^ 0.7 is expected 
(see Table 7 of Bell et al. I2003> . This would imply a stellar 
mass of A/host ~ 7 x lO^M©. The M/L-colour relations de- 
rived by Bell et al. f20031 are based on a large sample of galax- 
ies with no distinction of the galaxy type, so they represent 
an average over a wide morphological range. However, for a 
given colour, low-mass galaxies tend to show lower M/L val- 
ues (Binney & Merri field L?9§J. so 7 x lO^Mg is a plausi- 
ble upper limit to the host galaxy stellar mass. Therefore, we 
conclude that the total stellar mass of our galaxy should be 
lO^AfQ < Mho«t < 7 X IO^Mq. 

4. Discussion 

The results above were obtained assuming the i?23 low metal- 
licity branch (see Sollerman et al.|2005l. To verify that our re- 
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suits are not critically dependent on this assumption we have 
also run models with Solar abundances. 

No remarkable differences were found in the output param- 
eters. This is in agreement with the extensive tests performed 
by Bolzonella et al. J2000> who reported that the metallicity is 
a secondary parameter for the SED fitting technique 

The most relevant impact of the metallicity is the effect on 
the stellar age. Solar metallicity templates show slightly red- 
der colours than the Z — 0.004 templates, therefore shorter 
time evolution is required to fit the photometric points, yield- 
ing younger stellar ages. Thus, the mean stellar population age 
inferred using Solar metallicity templates was ^ 43 ± 7 Myr, 
a factor three lower than the one inferred with the Z = 0.004 
templates. 

The metallicity showed a second order impact on A^. Given 
that the Solar metallicity templates are redder, they provide 
lower Av values. Thus, assuming Solar metallicity we ob- 
tain an extinction of = 0.4 (instead of A^ = 0.6 for 
Z = 0.004), yielding an UV unextincted SFR=0.37 Mgyr^i 
and SSFR=23.1 Mq yr'^ (L/L*)"^. The rest of the parame- 
ters (photometric redshift, galaxy type, stellar mass) are basi- 
cally insensitive to the metallicity of the synthetic templates. 

4. 1. The GRB 030329 host SSFR relative to 
galaxies at a similar redshift 

We performed a comparative analysis of the GRB 030329 
host SSFR with respect to the SSFRs obtained for a galaxy 
sample based on the SDSS. In order to minimize any impact 
that the SSFR-redshift evolution (Cowie et al. 119961 Juneau 
et al. nM5[ might have on our analysis, only galaxies with 
0.160< z <0.175 were considered. Following the process ex- 
plained in C04A, the unextincted SSFR of the SDSS galaxies 
were calculated using HyperZ. We then selected only the galax- 
ies that showed acceptable SED fits (x^/d.o.f < 1.5). With 
the aim of estimating the probability errors, the process was 
repeated for the two IMFs given by S55 and CH03. 

Assuming a S55 IMF, among the 12333 SSFRs calculated 
for the SDSS sample, only 753 showed SSFRs larger than 
34 Mq yr~^{L/L*)~^, implying only 6.1% of the sample. The 
results did not change qualitatively when the CH03 IMF was 
used (6.9%). Hence we estimate that the GRB 030329 host 
SSFR is higher than ^^93. 5% of the galaxies at an equivalent 
redshift. 

This result is also supported by an independent (but more 
reduced) comparison sample based on 1067H-161 1 HDF galax- 
ies (adding the catalogues by Fernandez-Soto et al. 119991 and 
Vanzella et al. 2001). For a S55 IMF, only one of the 89 
HDF galaxies in the 0.05 < z < 0.25 redshift slot showed a 
SSFR larger than the one of the GRB 030329 host (98.9% per- 
centile). If a redshift 0< z < 0.25 window is assumed (con- 
taining 143 galaxies), then the GRB 030329 host SSFR is over 
97.9% of the sample. Again the conclusions remain basically 
invariant for a CH03 IMF (in a 0< z <0.25 redshift win- 
dow only 5 among 136 galaxies displayed SSFR > 34 A/0 
yr-i(L/L*)-i; 96.3% percentile). 



4.2. The GRB 030329 host in the framework of 
previous compilations 

The SFR (corrected for intrinsic extinction) of the GRB 030329 
host is the lowest among the C04A sample. However, the SSFR 
is the highest, given that the Mg of our host is 1.6 mag- 
nitudes dimmer than the faintest C04A sample galaxy. The 
GRB 030329 host falls into the selection criteria defined by 
C04A (R < 25.3, known redshift, detection in 5 or more filters, 
not very complex morphology). This allowed us to incorporate 
our host to the C04A 10 host galaxy sample and repeat part of 
the C04A statistical analysis. 

In order to construct a reference sample with similar spatial 
distribution to the (IOh-I) hosts, we calculated the SSFR only 
for the HDF galaxies with redshifts consistent with our 1 1 host 
galaxies (0.17 < z < 2.1). Selecting only HDF galaxies with 
acceptable SED fits (x^/d.o.f <1.5) reduced the HDF cata- 
logues to 1917 for a S55 IMF, and 1880 for a CH03 IMF 

The SSFR distribution of our 1 1 hosts and the HDF refer- 
ence samples were statistically compared using the two-sample 
Kolmogorov-Smirnov test (Peacock l 19831 see C04A for more 
details). The probability that the two samples share the same 
parent SSFR distribution ranges fromp = 1.2 x 10^'' (S55) to 
p = 5.5 X 10"^ (CH03) depending on the assumed IMF. This 
calculation reduces the C04A p values by an order of mag- 
nitude and strengthens the previous C04A results, who con- 
cluded that GRB hosts show on average higher SSFRs than 
HDF galaxies at similar redshifts. 

According to our synthetic SED fitting a iC'-band magni- 
tude of K' ^^20. 5 is predicted, implying R~K ^2.1 (adopting 
K = K' and correcting for Galactic reddening). This colour 
is consistent with the colours of the host galaxy sample by Le 
Floc'h et al. <2003t . The if -band absolute magnitude estimated 
for the host (Mk = —19.1) is in agreement with the general 
results by Le Floc'h et al. J2003> who noted that GRB hosts 
are significantly subluminous in the if-band (compare to their 
Fig. 5). 

The estimated R — K colour and the low redshift of the 
GRB 030329 host galaxy are consistent with the conclusions 
by Berger et al. ( 2003 ) and Tanvir et al. ( 2004 1, who found that 
the detected GRB hosts tend to be bluer and at lower redshifts 
than the submm-selected galaxies. In particular, it is illustrative 
of the position of the GRB 030329 host in the Fig. 6 colour- 
redshift diagram by Berger et al. ( 2003 ). 

On the other hand, the estimated host galaxy stellar mass is 
lower than the ones calculated for the host sample of Chary et 
al. J2002> . may be with the exception of the GRB 970508 host 
galaxy (Sokolov et aL ECRin . 

5. Conclusions 

VLT spectroscopic observations show that the GRB 030329 
host is a likely low metallicity galaxy, with Z — 0.004 
{Z ~ Zq/^). From the SED fitting we infer that the host of 
GRB 030329 is most probably a starburst galaxy. This agrees 
with the C04A results, who found (based on an independent 
host sample) that GRB hosts correspond to Stb galaxies in 
90% of the cases. We derived a dominant stellar population 
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Fig. 4. The unextincted SSFR of the C04A sample (filled dia- 
monds) plus the GRB 030329 host galaxy (triangle) as a func- 
tion of redshift. The crosses show the SSFR of the HDF when a 
S55 IMF is assumed. Note the high SSFR of the GRB 030329 
host in comparison to the HDF galaxies at a similar redshift. 



age of 150 ± 80 years, in agreement with the stellar age distri- 
bution of the C04A host galaxy sample (Mean age = 95 ± 74 
Myr). The low absolute magnitude of the host, Mb ^ —16.5 
(L ~ 0.Q16L*), confirms its subluminous nature. 

The SED fitting yields a consistent extinction value of 
Ay = 0.6 using either synthetic or empirical templates. This 
extinction is consistent with the mean extinction derived by 
C04A (median = 0.261^ Jg). The extinction derived from 
the Balmer emission lines is also consistent with the one de- 
rived from the SED. 

We determine a consistent unextincted SFR^ O.6M0 yr~^ 
based on three independent diagnostic techniques; the UV con- 
tinuum, the [On], and the Balmer emission lines. This con- 
sistency supports the inferred A^ — 0.6 value. The associated 
unextincted SSFR (- 34Afo yr-'^{L/L*)-^) is higher than 
93.5% of the SDSS galaxies at z 0.17. Moreover, a sta- 
tistical analysis carried out enhancing the C04A sample with 
the GRB 0303029 host confirms that GRB hosts shows higher 
SSFRs than HDF galaxies at similar redshifts. All these find- 
ings are consistent with the host galaxy being an active star 
forming galaxy. 

We note the presence of a close galaxy with a similar red- 
shift to the GRB host, which may be indicative of a possi- 
ble clustering at z ^0.17. Multi-object spectroscopy of the 
GRB 030329 field may clarify this open issue. 
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